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■ ABSTRACT 

o 

o . 

■ We report the abundances of neutron-capture elements in eight carbon-rich, metal- 
poor (—2.7 <[Fe/H]< —1.9) stars observed with the Subaru Telescope High Dispersion 

^ . Spectrograph. The derived abundance patterns indicate that the neutron-capture el- 

ements in these objects primarily originated from s-process nucleosynthesis, although 
the [Ba/Eu] abundance ratios in some objects are lower than that of the solar-system 
^ _ s-process component. The present analysis has yielded the Pb abundances for seven ob- 

' jects, as well as an upper limit for one object, from use of the Pb I 4057 A and 3683 A 

■ lines. The values of [Pb/Ba] in these objects cover a wide range, between —0.3 and 

■ +1.2. Theoretical studies of s-process nucleosynthesis at low metallicity are required to 
. explain this large dispersion of the [Pb/Ba] values. Variations in radial velocity have 

Q . been found for two of the eight objects, suggesting that, at least in these instances, the 

observed excess of s-process elements is due to the transfer of material across a binary 
Qh| system including an AGB star. Comparisons with predictions of AGB nucleosynthesis 

O ' models are discussed. 
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1. Introduction 

Lead (Pb) isotopes form, along with those of bismuth, the group of the heaviest stable nuclei. 
Understanding of the physical processes responsible for the synthesis of these nuclei, as well as their 
enrichment history in the Galaxy, is important in the study of stellar nucleosynthesis and Galactic 
chemical evolution, as well as for precision cosmochronology, since these nuclei are the expected 
decay products of Th and U (Goriely & Arnould 2001; Schatz et al. 2002). The majority of Pb 
nuclei in the solar system are believed to have originated in the s-process. However, the so-called 
main- component of the s-process in classical models cannot explain the solar system abundances 
of Pb, hence another component (the so-called strong component) was introduced (e.g., Kappeler, 
Beer, & Wisshak 1989). Recent theoretical work on nucleosynthesis in asymptotic giant branch 
(AGB) stars, based on the idea of an s-process that occurs in the radiative zone during the inter- 
pulsc phase, have predicted large production of Pb and Bi (Z = 82 and 83, respectively) compared 
with lighter s-proccss elements such as Sr (Z = 38) and Ba (Z = 56) in low mctallicity AGB stars 
(Gallino et al. 1998; Goriely & Mowlavi 2000; Busso et al. 2001). These works suggested that 
the Pb production in low-metallicity AGB stars corresponds to the strong component assumed in 
the classical models. Studies of Pb enrichment in the Galaxy have also started using the yields 
predicted by these AGB models (Travaglio et al. 2001). 

Recent abundance studies of s-process-element-enhanced, metal-poor stars, based on high- 
resolution spectroscopy, have measured the abundances of Pb and other neutron-capture elements 
in a number of stars, and made it possible to examine the models of AGB nucleosynthesis in greater 
detail. For instance, Aoki et al. (2000) determined the abundances of 16 neutron-capture elements, 
including Pb, for the s-process-element-rich, very metal-poor ([Fe/H]=— 2.7) star LP 625-44®. The 
Pb abundance of this object is, however, much lower than the prediction by the "standard" model 
of Gallino et al. (1998), when the abundance pattern is normalized by the abundance of lighter 
s-process nuclei (e.g., Ba). A similar result was derived for another metal-poor star, LP 706-7, by 
Aoki et al. (2001). Their abundances could, however, be reconciled with the models of Gallino 
et al. if the extent of the ^^C pocket is reduced by a factor of ~24 (Ryan et al. 2001). On 
the other hand, van Eck et al. (2001) studied another three s-process-element-enhanced, metal- 
poor objects, and found very large Pb enhancements in these objects, compared to those of the 
lighter elements. They concluded that the results can be explained well by the model of AGB 
nucleosynthesis presented by Goriely & Mowlavi (2000). 

Thus, the abundance ratios of the Pb, and the neutron-capture elements at the second peak 
of the s-process (Z ~ 56), produced by metal-deficient AGB stars seem to show a large dispersion, 
or might even be classified into two distinct groups. In order to understand the nature of the 
s-process nucleosynthesis at low metallicity, we have embarked on an extensive set of studies of the 
neutron-capture elements for metal-poor stars with excesses of s-process elements. In this Paper 
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we report the abundances of neutron-capture elements, including Pb, for eight metal-poor stars. 

2. Observations and Measurements 

We selected our program sample from the list of metal-poor stars provided by the HK survey 

of Beers and collaborators (Beers, Preston, & Shectman 1992; Beers 1999; Norris, Ryan, & Beers 
1999), taking account of the metallicity and apparent level of carbon enhancement (as indicated by 
an index, GP, that measures the strength of the CH G band at 4300 A ) . We have observed more 
than 10 stars, and selected the 7 stars that clearly exhibit excesses of the s-process elements. For 
purposes of comparison, we have also observed the s-process-element-rich star HD 196944 studied 
by van Eck et al. (2001). 

Observations were carried out with the High Dispersion Spectrograph (HDS) of the 8.2 m 
Subaru Telescope (Noguchi et al. 2002). Program stars and observational details are provided in 
Table 1. The spectra were taken with a resolving power R = 50, 000, except for HD 196944, for 
which a R = 90, 000 spectrum was obtained. Our spectra cover the wavelength range 3550-5250A. 
Signal-to-Noise (S/N) ratios per 0.012 A pixel at 4000 A are given in the table as well. Data 
reduction was performed in the standard way within the IRAF^ environment, following Aoki et al. 
(2002b). 

Equivalent widths for isolated lines were measured by fitting Gaussian profiles to the absorption 
lines. The line list and equivalent widths measured are presented in Table 2. 

The radial velocity (Vr) measured for our spectra are also presented in Table 1. Our separate 
measurement of the radial velocity for CS 29526-110 {V, = 201.70 ± 0.26km s"^ at JD=2452152.7) 
using WHT/UES shows a variation of the radial velocity, suggesting binarity for this object. The 
radial velocities measured for our Subaru/HDS spectra of CS 22880-027 and CS 22898-027 agree 
with the averages of the results by Preston & Sneden (2001) to within 1 km s~^. 

3. Abundance Analysis 

A standard analysis, using the model atmospheres of Kurucz (1993) and the equivalent widths 
measured above, was performed for lines of Fe I, Fe II, and for most of the neutron-capture elements. 
A spectrum synthesis technique was applied to the molecular bands of CH, C2, and CN in order 
to determine the carbon and nitrogen abundances, which are required for estimates of the effective 
temperature. The line broadening due to macro-turbulence in the stellar atmospheres, and from 
the spectrograph itself, was estimated from the clean Fe lines for each object, and then applied to 
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the spectrum synthesis. 

The effective temperatures given in Table 3 were estimated from broadband colors using pre- 
viously derived temperature scales. For most of our program stars, we applied the temperature 
scale based on the de-reddened (B — V)o color for carbon-rich, metal-poor stars produced by Aoki 
et al. (2002a). The B — V color of CS 30301-015 is uncertain, hence we assume a larger uncer- 
tainty in effective temperature in the error estimation procedure described below. The effective 
temperature of CS 22942-019 was re-determined after its carbon and nitrogen abundances were 
estimated, because abundances for these two elements assumed in the first determination of the 
effective temperature were found to be too low, and the effect of molecular absorption was not 
included correctly. JHK photometry data from the interim release of the 2MASS Point Source 
Catalog (Skrutskie et al. 1997) are also available for three stars (CS 29526-110, CS 22898-027 and 
CS 31062-050). After correcting for reddening, the effective temperature scale based on [V — K)q 
by Alonso, Arribas, & Martinez-Roger (1996) was applied, and yielded consistent results with those 
obtained from {B — V)o colors. The difference between the two effective temperatures is rather 
large in CS 29526-110: the effective temperature estimated from {V — K){y is higher by 200 K than 
that from [B — V)q. We assume a larger uncertainty of T^s for this object in the error estimates 
below. 

Surface gravities for our program stars were determined from the ionization balance between 
Fe I and Fe II, the metallicities were estimated from the abundance analysis for these lines, and 
the micro-turbulence was determined from the Fe I lines by demanding no dependence of the 
derived abundance on equivalent widths. The results are provided in Table 3. While the effective 
temperature of CS 29526-110 is the highest in our sample, and suggests that this star might be a 
main-sequence turn-off star, its gravity (log g = 3.2) is not particularly high. This may be due to 
non-LTE effects in the analysis of the Fe lines, which were neglected in the present work. For the 
purpose of error estimation, we therefore assume a larger uncertainty in the derived gravity for this 
star. 

The carbon abundances were determined using the C2 Swan bands for CS 31062-050, CS 22942- 
019, and CS 30301-015, stars for which the CH band at 4323 A is too strong for abundance analysis. 
The line list produced by Aoki k, Tsuji (1997) for C2 Swan bands were applied. For the other five 
stars, the CH band was used for carbon abundance determination. The line positions of CH A — X 
bands were calculated using the molecular constants derived by Zachwieja (1995) and the oscillator 
strengths of Brown (1987). Aoki et al. (2002a) showed that the carbon abundance derived from 
the C2 band is by 0.2 dex higher than that from the CH band for LP 625-44, indicating that we 
should include this additional uncertainty of about 0.2 dex for the carbon abundances. The CN 
3883 A feature was used for estimates of nitrogen abundances. The line positions of the CN violet 
system were calculated using the molecular constants presented by Prasad k, Bernath (1992) and 
the oscillator strengths of Bauschlichcr, Langhoff & Taylor (1988). Unfortunately, the nitrogen 
abundances of CS 29526-110, CS 22880-074, and CS 30301-015 are unreliable, because the CN 
band is either too weak or too strong for abundance determination. The results of the carbon and 
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nitrogen abundances are also given in Table 3. The oxygen abundances are assumed to be [0/Fe] = 
0.5 in the analysis, taking account of the well-known overabundances of oxygen found in metal-poor 
stars. The derived carbon and nitrogen abundances are not sensitive to the oxygen abundances 
assumed in the analysis for warm (Teff > 5000 K) stars, because the fraction of carbon bound in 
the CO molecule is small. The effect of the assumption for the oxygen abundance most severely 
appears in CS 30301-015 (Teg = 4750 K); we tested the range 0.0 <[0/Fe]< 1.5 for this star, but 
found that the effect on abundance determination for carbon and nitrogen is smaller than 0.1 dex. 

The Pb I 4057 A line was detected for seven stars in our sample. Figure 1 shows examples of 
the observed spectra, and the results of synthetic fits, for HD 196944 and CS 31062-050. In this 
analysis we adopted the line data for the hyperfine splitting and isotope shifts determined by Simons 
et al. (1989), as given in Table 4. We found that the fit between the observed spectra and the 
calculated ones is better, when the absolute wavelength of the Pb lines are slightly shifted. For this 
reason, we shifted the wavelength of these lines by 0.005 A, which is as large as the uncertainty in 
the determination of the absolute line position by Simons et al. (1989), in our spectrum synthesis. 
The CH 4057 A line blending with the Pb I line was calculated for the carbon abimdance which 
best reproduced other lines of the same CH band. Therefore, the errors in carbon abundances due 
to the uncertainties of oscillator strengths of the band system do not directly affect the spectrum 
synthesis for the Pb I 4057 A line. 

The effects of hyper-fine splitting and isotope shifts were included, assuming the solar system 
isotope ratio for Pb (see Aoki et al. 2001). However, this assumption may not be appropriate for 
the very metal-poor stars studied here. For instance, a large production of ^^^Pb was predicted for 
the s-process at low metallicity (e.g., Gallino et al. 1998). In order to estimate the uncertainty of 
the Pb abundance due to the effect of the isotope ratio assumed, we also tried the analysis with a 
single line approximation for the Pb I 4057 A line, which simulates the case that ^"^^Pb is dominant 
over the other isotopes. The synthetic spectra calculated with the single line approximation are 
shown by dotted lines in Figure 1. Because of the larger effect of saturation of the absorption, 
the derived Pb abundance from the single line approximation is larger than that from the analysis 
including isotope shifts with the assumption of the solar system isotope ratio. The difference of the 
resulting abundance is about 0.2 dex for CS 31062-050, in which the Pb I 4057 A line is strongest 
among our objects. The effect of the isotope ratios assumed in the analysis for other objects is 
about 0.1 dex or smaller. We include these uncertainties in the error estimates below. 

We found that the abundances of the lighter neutron-capture elements (e.g., Ba), as well as the 
stellar parameters, of CS 31062-050 are quite similar to those of the well-studied object LP 625-44. 
However, the Pb absorption line in CS 31062-050 is significantly stronger than that of LP 625-44 
(see Figure 6 of Aoki et al. 2002b). This clearly indicates a large dispersion in the ratio of the 
abundances between the lighter s-process elements and Pb. 

The Pb I 3683 A line was also detected in five of the above seven stars. Though the quality of 
the spectra at 3683 A is worse than that at 4057 A, we found that the abundances derived from the 
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3683 A line agree well with those obtained from the 4057 A line within the errors. We note that 

the abundances based on the 4057 A line were adopted as the final result given in Table 5; those 
based on the 3683 A line were only used for confirmation. An upper limit on the Pb abundance 
was obtained for CS 22942-019, in which the above Pb lines were not detected. 

The abundances of most of the other neutron-capture elements were determined by standard 
analysis based on the measured equivalent widths, while the abundances of Ba and Eu were deter- 
mined from spectrum synthesis. The effects of hyper-fine splitting and isotope shifts were included 
in the analysis of Ba, La, Nd, and Eu lines using the line lists of McWilliam (1998) and Aoki et 
al. (2001). 

In order to estimate the random error in the analysis, we first calculated the dispersion of the 
abundances for individual Fe I lines in each star. In addition, we assumed the random error in 

g/- values to be 0.1 dex, and added it in quadrature to the above dispersion. The random error 
was evaluated by dividing this estimate by n^/^ (n is the number of lines used in the analysis) 
for each neutron-capture element. For Pb abundances, we also included the uncertainties due to 
the isotope ratio assumed in the spectrum synthesis (0.2 dex for CS 31062-050, and 0.1 dex for 
other stars) . The errors in the abundance determination from the uncertainties of the atmospheric 
parameters were evaluated for a{Tes) = lOOK, ^(log^) = 0.3, and a{v) = 0.5 km for HD 196944 
and CS 22898-027. These results were then applied to four of the remaining stars. We assumed 
a{Tcs) = 200 K for CS 29526-110 and CS 30301-015, and cr{logg) = 0.5 for CS 29526-110, taking 
the uncertainties noted above into consideration. We finally derived the total uncertainty by adding 
in quadrature the individual errors, and list them in Table 5. 

We also tried to estimate the effective temperatures for seven stars in our sample based on 
the Balmer line strengths measured from the medium-resolution spectra obtained during the HK 
survey follow-up. These effective temperatures are systematically lower, by about 300 K, than 
those derived above from the broadband colors. The effects of the changes of effective temperature 
by — 300K on the resulting [Fe/H] and [Pb/Ba], which will be discussed below, arc about —0.20 and 
+0.25 dex, respectively. Though these effects are important in the determination of the abundances, 
and will be investigated more thoroughly in our future work, the systematic shift of the abundance 
ratios does not change the following discussion on the dispersion in the values of [Pb/Ba] found in 
our sample. 

The abundances of HD 196944 derived in this work generally agree with those by Zacs, Nissen 
Sz Schuster (1998) for iron and carbon, and van Eck et al. (2001) for other elements, including Pb, 
within 0.2 dex, except for La (our result is 0.38 dex higher than that of van Eck et al. (2001)). The 
agreement between the abundances of HD 196944 determined by van Eck et al. (2001) and those 
obtained in the present analysis (following the previous work for LP 625-44) indicates that the 
significant difference in the abundance ratios of Pb, as compared to the lighter s-process elements, 
in the objects analyzed by van Eck et al. (2001) and those by Aoki et al. (2001) is not due to 
differences in the respective analysis procedures. 
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4. Discussion and Concluding Remarks 

The resulting abundances for two of the stars in our sample are shown in Figure 2, along with 
the scaled abundance patterns of solar system material, the main s-process component, and the 
(inferred) r-process pattern (Arlandini et al. 1999). These comparisons are useful to distinguish 
the origins of the neutron-capture elements detected in our objects. The abundance patterns of 
elements with 56 < Z < 63 in our stars agree with the s-process pattern much better than with 
the r-process pattern. This result implies that, as expected, the neutron-capture elements in our 
stars principally originate in the s-process. The values of [Ba/Eu] for four of the stars in our 
sample (CS 29526-110, CS 22898-027, CS 31062-012, and CS 31062-050) are significantly lower 
( [Ba/Eu] =0.36~0.47) than seen in the main s-process component ([Ba/Eu]=-|-1.15; Arlandini et al. 
1999). We suggest, however, that the derived abundance ratios of [Ba/Eu] have been produced by 
an s-process which produces different abundance ratios from that of the main s-process component. 
Indeed, Goriely & Mowlavi (2000) predicted lower values of [Ba/Eu] (~ 0.4) for yields of metal- 
deficient AGB stars, though the [Ba/Eu] value of our objects is not correlated with metallicity. To 
explain these low values of [Ba/Eu] by the mixture of the abundance ratios of the r- and s-process 
components in the solar system, we must assume that about 80-90% of Eu nuclei were produced 
by the r-process. If this were true, it follows that these four stars show very large excesses of their 
r-process elements ([Eu/Fe] = 1.5~1.8), similar to CS 31082-001, an extreme r-process-element- 
enhanced star (Cayrel et al. 2001; Hill et al. 2002). However, recent high-resolution studies of 
metal-poor giants indicate that stars with such large excesses of r-process elements are quite rare. 
An ongoing survey by Christlieb et al. (2002, in preparation) confirms the original suggestion by 
Beers (private communication) that roughly 3% of giants with [Fe/H] < —2.5 exhibit [r/Fe]> -1-1.0. 
Hence, it would be difficult to appeal to this explanation to account for the excesses of Eu in the 
four stars in our sample. We treat our stars as having the abundance patterns produced by the 
s-process for neutron-capture elements in this paper, because the elements discussed below (Ba and 
Pb) should principally originate from the s-process, even though r-process contamination may have 
contributed to a portion of the observed Eu excesses. 

Figure 3 shows the ratio [Pb/Ba], representing the ratio of abundances between elements 
located at the second and third s-process peaks, as a function of metallicity ([Fe/H]). There exists 
a large scatter in [Pb/Ba] for these stars, larger than can be accounted for by the errors in the 
abundance analysis. It is not apparent that our stars can be readily classified into separate groups 
on the basis of the present data. There may be evidence for an decreasing trend of [Pb/Ba] with 
decreasing [Fe/H] (the correlation estimated from objects in Figure 3, excluding CS 22942-019, 
is represented as [Pb/Ba]= 2.62 -|- 0.88[Fe/H] with the standard error of the slope of 0.47, i.e., 
the possible slope is significant at the 1.9 sigma level). Additional data will presumably help to 
clarify these questions. Additional studies of s-process nucleosynthesis at low metallicity should be 
carried out to explore the possible reasons for this large dispersion, and the possible dependence 
on metallicity, of the Pb/Ba ratios. 

Comparison of these results with the predictions of AGB nucleosynthesis models of Gallino et 
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al. (1998) and Busso et al. (1999) provides a constraint on their ^^C pocket models®. Figure 1 
of Ryan et al. (2001) shows the metallicity dependence of [Pb/Fe] and [Ba/Fe] predicted by the 
model for 1.5 Mq AGB stars, compared with the abundance ratios observed in LP 625-44. The 
value of [Pb/Ba] is sensitive to the adopted ^^C profile, represented by the normalization factor of 
the -"^^C pocket in the standard model of Busso et al. (1999). The [Pb/Ba] values found in our stars 
with —2.7 <[Fe/H]< —1.9 distribute from —0.3 to +1.2 (Figure 3). This range of the abundance 
ratios is not explained by the standard model, which predicts [Pb/Ba] ~ +2, but can be explained 
by models where the amount of ^^C produced as the neutron source for the s-process is lowered by 
factors of 6-24 relative to the standard model. 

We would also like to mention the model of the s-process in AGB stars recently proposed by 
Iwamoto et al. (2002). Their model for AGB stars with [Fe/H] = -2.7 and M = 2Mq shows 
a possible production of ^^C, a candidate for the neutron source of the s-process, as the result of 
proton mixing into the hot He shell at the second thermal pulse (see also Fujimoto et al. 2000). 
Since the value of neutron exposure predicted by their model is rather high (r ~ 1 mb~^), one may 
expect high Pb abundances relative to those of lighter neutron-capture elements. However, the small 
number of the exposures (only one in their model) is preferable to account for the low Pb abundances 
that we find in some stars in our sample (Aoki ct al. 2001). This mechanism for producing the 
required free neutrons is expected to occur only in very metal-poor ([Fe/H]< —2.5) AGB stars, 
and may explain the abundances of some objects with lower metallicity and high [Pb/Ba] values. 
The comparison with this model is, however, still quite speculative, and further theoretical study 
for the possibility of flash-driven proton mixing and its contribution to the s-process is required. 

Finally, we point out that clear temporal variations of radial velocity, which directly implies the 
binarity of the object, have thus far been detected for only three of the carbon-enhanced, metal-poor 
stars we have been studying in recent years: CS 22942-019 (Preston &: Sneden 2001), CS 29526- 
110 (section 2), and LP 625-44 (Aoki et al. 2000). Preston & Sneden (2001) reported that none 
of the three carbon-rich, metal-poor subgiants in their sample, including the stars CS 22898-027 
and CS 22880-074 studied here, exhibited radial velocity variations exceeding 0.5 kms~^ over an 8 
year period. Furthermore, no radial velocity variation has been found for the subgiant (or dwarf) 
LP 706-7 (Norris, Ryan, & Beers 1997). These results suggest that the enhancement of neutron- 
capture elements in these subgiants may not be explained as a result of the transfer of material 
rich in s-process elements across a binary system including an AGB star. For comparisons with 
the predictions of AGB models, confirmation of binarity based on the radial velocity monitoring 
is strongly desired. It is interesting to note the large dispersion of [Pb/Ba] values of the three 
objects (CS 22942-019, CS 29526-110, and LP 625-44) for which the variation of radial velocity 
have been observed. This suggests that the Pb/Ba ratios produced by s-process nucleosynthesis in 



Uncertainties in the amount of C in the pocket reflects the unknown amount of mixing of protons from the 
H-rich envelope down into the C-rich zone. Recent model calculations by Cristallo et al. (2001; see also Goriely & 
Mowlavi 2000) treat the proton mixing as an input paxameter, rather than the amount of ^^C in the "pocket". The 
efficiency of the ^^C as a neutron source for the s-process is also affected by the profile of neutron poisons like ^'*N. 
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metal-deficient AGB stars show a large scatter, although the sample is still too small to be confident 
in this regard. 

The present work clearly shows a variety in the abundance ratios of neutron-capture elements 

for carbon-rich, metal-poor stars, and underscores the importance of further studies for individual 
objects, including long-term radial velocity monitoring, for development of a better understanding 
of heavy element production by the s-process operating at low metallicity. 

The authors are pleased to acknowledge valuable discussions with Drs. R. Gallino and M. 
Busso on the variation of s-process nucleosynthesis in their Population II AGB models. We would 
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Table 1. PROGRAM STARS AND OBSERVATIONS 



Star 




V 


B-V 


E{B - V) 


GP Exp.^ 


S/N 




Obs. Date (JD) 


CS22942- 


-019 


12.71 


0.86 


0.02 


7.60 60 (2) 


61 


-232.88±0.40 


20 Aug., 2000 (2451778) 


CS29526- 


-110 


13.37 


0.35 


0.00 


1.94 90 (3) 


62 


-h208.35±0.76 


28 Jan., 2001 (2451938) 


CS22898- 


-027 


12.76 


0.50 


0.03 


4.71 60 (3) 


49 


-47.90±0.47 


22 July, 2001 (2452113) 


HD 196944 


8.41 


0.61 


0.00 


- 20 (1) 


150 


-174.76±0.36 


25 July, 2001 (2452116) 


CS31062- 


-012 


12.10 


0.46 


0.00 


3.67 25 (1) 


57 


+80.76±0.37 


25 July, 2001 (2452117) 


CS31062- 


-050 


13.04 


0.64 


0.02 


6.41 50 (2) 


47 


+7.75±0.67 


25 July, 2001 (2452117) 


CS30301- 


-015 




0.99: 




7.88 50 (2) 


53 


+86.46±0.43 


25 July, 2001 (2452116) 


CS22880- 


-074 


13.27 


0.57 


0.04 


5.71 60 (2) 


42 


+59.57±0.48 


25 July, 2001 (2452116) 



^ Exposure time (minute) and number of exposures. 
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Table 2. EQUIVALENT WIDTHS 



Wavelength 

(A) 


LEP 

(eV) 


logy/ 


29526 
-110 


22898 
-027 


31062 
-012 


22880 
-074 


31062 
-050 


HD 196944 


22942 
-019 


30301 
-015 


Sr II 
























4077.71 





00 


0.15 


102.3 


115.4 


76.6 


127.6 


146.9 


179.7 


488:t 


176:t 


4161.80 


2 


94 


-0.50 












16.1 






4215.52 





00 


-0.17 


90.3 


97.8 


56.2 


104.6 


127.2 


184.7 


364:t 




Y II 
























3600.73 





18 


0.28 




48.6 








95.4* 






3611.04 





13 


0.01 








34.2 




86.0* 






3628.70 





13 


-0.7P 












58.7* 






3788.69 





10 


-0.06 












96.8 






3950.35 





10 


-0.49 




25.0 


7.6 


30.2 




78.7 


101.0 




4682.32 





41 


-1.5P 












22.8 






4883.68 


1 


08 


0.07 












56.2 


103.1 


61.2 


5087.42 


1 


08 


-0.16 












42.7 




55.2 


5200.41 





99 


-0.57 












36.6 




33.0 


5205.72 


1 


03 


-0.34 












35.7 




41.8 


Zr II 
























3573.08 





32 


-1.04 












33.5 






3630.02 





35 


-1.11 












25.1 






3998.97 





56 


-0.67 












51.2 


75.4 




4029.68 





71 


-0.76 




8.0 






23.6 


41.9 






4050.33 





71 


-1.00 














59.2 




4149.20 


n 
u 


OU 


-0.03 


28.7 
















4208.98 





71 


-0.46 




16.9 






31.9 


47.0 






4258.05 





56 


-1.13 












29.9 






4443.00 


1 


49 


-0.33 














56.4 




Ba II 
























3891.78 


2 


51 


0.28'' 




49.5 








34.0 


50.7t 




4130.65 


2 


72 


0.56^= 


50.2 


64.3 


23.1 


31.0t 


85.2 


49.1 


69.4t 


81.7 


4166.00 


2 


72 


-0.41 


12.7t 


15.1 






27.7 


11.5 


22.lt 




La II 
























3790.83 





13 


0.03 










70.7 








3949.10 





40 


0.49 


69.0 


107.9 


34.4 


50.4 


129.1* 


95.1 






3988.52 





40 


0.21 


44.3 






33.3 


120.9* 


67.9 




70.3 


3995.75 





17 


-0.06 


34.3 


59.5 


18.6 




89.4* 


54.2 


60.4 


73.1 
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Table 2 — Continued 



4086.72 


0.00 


-0.07 




54.3 




33.0 


73.7 




61.1 




4123.23 


0.32 


0.13 


40.3 




19.1 




74.2 


50.9 






4322.51 


0.17 


-0.93 




20.7 






51.7 


15.8 






4333.74 


0.17 


-0.06 


54.1 
















4429.90 


0.23 


-0.35 


36.5 










47.6 






4662.51 


0.00 


-1.24 




14.4 












32.7 


Ce II 






















3853.16 


0.00 


-0.36"^ 




21.5 














3909.31 


0.45 


-0.52 










19.3 


9.0 






3940.34 


0.32 


-0.20*^ 




16.6 






25.8 


24.7 




30.1 


3942.75 


0.86 


0.41^ 




24.8 








23.9 




33.7 


3960.91 


0.32 


-0.20 












15.4 






3980.88 


0.71 


0.03 












12.0 






3984.68 


0.96 


0.41 












15.2 






3992.39 


0.45 


-0.13*^ 




10.4 














3999.24 


0.29 


0.32"^ 






8.5 


17.5 


45.6 


32.0 


39.8 


53.0 


4014.90 


0.53 


0.05^^ 




12.3 






34.8 








4046.34 


0.55 


-0.11^ 




8.9 








15.8 






4053.51 


0.00 


-0.60^^ 










26.1 


17.7 


23.8 




4068.84 


0.70 


-0.13*^ 












16.6 






4072.92 


0.33 


-0.71 




11.0 








14.7 






4073.48 


0.48 


0.28"^ 




27.2 




19.9 


35.6 


31.7 






4118.14 


0.70 


O.lO'i 










37.9 


23.9 






4120.83 


0.32 


-0.51^^ 












19.4 






4124.79 


0.68 


-0.10*^ 










30.7 








4127.37 


0.68 


0.19^^ 










37.1 








4137.65 


0.52 


0.32^1 




34.8 


11.3 




49.9 


39.2 




59.3 


4145.00 


0.70 


0.03^1 












14.6 


20.5 




4165.61 


0.91 


0.51^^ 




22.8 






38.4 


28.8 






4185.33 


0.42 


-0.56 












10.6 






4193.87 


0.55 


-0.40'^ 










18.5 


14.0 






4222.60 


0.12 


-0.21*^ 




28.9 


8.8 




44.6 


33.5 




58.1 


4349.79 


0.53 


-0.19^1 










25.6 








4418.78 


0.86 


0.26^^ 


12.7 


23.9 








30.4 






4449.34 


0.61 


-O.ll'i 




21.6 








25.2 






4460.21 


0.48 


0.25^^ 




35.1 








40.0 




77.8 
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Table 2 — Continued 



4471.24 


0.70 


0.23^ 






10.6 




44.9 


33.4 






4483.90 


0.86 


0.09 














47.1 




4486.91 


0.29 


-0.39^1 




18.6 






39.7 


31.7 


35.5 


46.1 


4551.30 


0.74 


-0.49 










19.2 






24.5 


4560.28 


0.91 


0.08^ 




17.4 






37.3 


19.9 




53.6 


4560.96 


0.68 


-0.47 










22.7 


9.2 




23.7 


4562.36 


0.48 


0.16^ 






9.9 


18.8 


57.5 


36.0 


49.0 


70.5 


4593.93 


0.70 


-O.OS'^ 




19.0 


11.9 




49.6 


19.9 




55.8 


4628.16 


0.52 


0.09^ 




18.6 


10.2 




48.5 


29.2 




57.4 


5187.45 


1.21 


0.16<i 
















14.2 


Nd II 






















3780.40 


0.47 


-0.27 












25.2 






3784.25 


0.38 


0.23 




34.1 


7.0 






30.7 


28.8 




3911.17 


0.47 


0.20 










52.2 


27.6 






3927.10 


0.18 


-0.52 










33.0 




18.2 




3941.51 


0.15 


0.06 






14.8 






36.5 






3979.49 


0.20 


-0.11 












29.5 


30.1 




3991.74 


0.00 


-0.40 




28.9 






54.8 


33.0 


50.8 




4012.70 


0.00 


-0.64 










36.3 


18.0 


20.8 




4020.87 


0.32 


-0.18 










47.1 


29.1 


47.1 




4021.34 


0.32 


0.24 






7.3 




43.2 


21.3 






4023.00 


0.20 


-0.19 




23.4 






50.8 


18.5 






4038.12 


0.18 


-0.78 










44.4 








4051.15 


0.38 


-0.21 




19.9 






34.5 








4059.96 


0.20 


-0.33 
















17.5 


4061.09 


0.32 


-0.34 












48.0 






4069.28 


0.06 


-0.33 




22.7 






39.1 


21.1 


16.0 


43.3 


4075.12 


0.20 


-0.40^ 


11.6 


35.6* 


22.4* 












4116.77 


0.06 


-1.65 










24.7 








4133.36 


0.32 


-0.34 












17.8 






4156.08 


0.18 


0.10^ 


28.4 


43.2* 














4160.57 


0.56 


-0.45 










25.6 








4177.32 


0.06 


-0.04'^ 




38.7* 




22.5 


61.7* 








4358.17 


0.32 


-0.21 




30.9 














4446.39 


0.20 


-0.63 




19.0 






44.9 


29.7 




37.8 


4462.99 


0.56 


-0.07 




35.4 


10.6 




56.9 


28.5 
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Table 2 — Continued 



4516 36 


0.32 


—0.75 










26.1 


8.4 






Sm Tl 






















SQ7Q 90 


54 


—0 1QS 












7 1 


1 9 6 

J- t/.U 


1 3 8 

±0.0 


4064 58 


0.33 


— 0.27s 










19.5 








4944 70 


98 


—0 73^ 










14 4 










98 


—0 97 




1 1 7 






97 fi 


1 5 1 






4424.34 


0.49 


0.07s 












11.9 




28.0 


4434.32 


0.38 


-0.26 




16.0 






33.8 


14.7 






4467.34 


0.66 


0.12' 




13.1 






30.4 


12.5 




27.2 


4499.48 


0.25 


-1.00*^ 










23.6 








4587 95 


0.49 


—0.23 










30.5 








4642.24 


0.38 


—0.52*^ 










33.1 








Rii IT 






















381 Q 67 


00 


4Qj 


55 


83 5t 


19 it 

J. C/. -L 




132.4t 


43 4t 


60 7t 




3Q07 1 


0.21 


20 










87 4t 


15 4t 






41 99 70 


00 


90j 


44 fit 


55 7t 


1 9 8t 




116 Ot 

± ±\J.\J 


99 4t 


40 5t 


47 it 

^ 1 . J- 


4905 05 


on 


n 1 9j 


7t 
00. i 


^0.0 


R 7t 
0. 1 


Q 8t 




91 1 1 

Zi J. . J. 


'^'^ nt 

0(J.U 




T^-rr TT 

Dy 11 






















3944.68 


0.00 


0.08 




28.3 




11.9 




32.7 


34.9 


48.1 


4077.96 


0.10 


-0.03 




31.0 






53.4 








Er II 






















3616.56 


0.00 


-0.33'^ 




27.9 




19.0 




23.4 






3786.84 


0.00 


-0.64 




27.4 














3830.48 


0.00 


-0.37'^ 




43.7 














3896.23 


0.05 


-0.24"^ 




36.5 








29.5 






Pb I 






















3683.46 


0.97 


-0.52 


30.lt 


24.2t 






45.3t 


24.lt 




26.2t 


4057.82 


1.32 


-0.20 


34.9t 


27.4t 


6.4t 


15.4t 


53.3t 


24.5t 




20.7t 



* Lines which were not used in the analysis, 
f Synthesized value calculated for the abundance derived by spectrum synthesis. 
References for ^'/-values. — ^Hannaford et al. (1982), ^Sneden et al. (1996), '^McWilliam 
(1998), Corliss and Bozman (1962) +0.23 dex, ^ Cowley and Corliss (1983), ^ Wiese and 
Martin (1980), § Corliss and Bozman (1962)+0.49 dex, ^ Biemont et al. (1989), ' Vogel et al. 
(1988), j Biemont et al. (1982), ^ Musiol and Labuz (1983), others: Aoki et al. (2002b) or Aoki 

et al. (2001) 



-17- 



Table 3. STELLAR PARAMETERS 



Star 


Teff(K)* 




t;tur(kms-^) 


[Fe/H] 


[C/Fe] 


[N/Fe; 


CS 29526-110 


6500^^ 


3.2 


1.6 


-2.38 


2.2 


1.4: 


CS 22898-027 


6250'''*' 


3.7 


1.5 


-2.25 


2.2 


0.9 


CS 31062-012 


6250^ 


4.5 


1.5 


-2.55 


2.1 


1.2 


CS 22880-074 


5850^ 


3.8 


1.4 


-1.93 


1.3 


-0.1: 


CS 31062-050 


5600^'^ 


3.0 


1.3 


-2.32 


2.0 


1.2 


HD 196944 


5250^ 


1.8 


1.7 


-2.25 


1.2 


1.3 


CS 22942-019 


SOOO'^ 


2.4 


2.1 


-2.64 


2.0 


0.3 


CS 30301-015 


4750^ 


0.8 


2.2 


-2.64 


1.6 


0.6: 



* The effective temperature is determined based on the foUowing color indices: a:{B — V)o and 

h:{V-K)o. 



Table 4. Hyperfine and isotope splitting of Pb 4057.8 A line 



Wavelength (A) fraction isotope 



4057.7524 


0.015 


207b 


4057.7895 


0.135 


207a 


4057.7895 


0.524 


208 


4057.8030 


0.236 


206 


4057.8150 


0.015 


204 


4057.8251 


0.075 


207c 
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Table 5. [FE/H] AND RELATIVE ABUNDANCE, [X/FE] 





CS 29526- 


-110 


CS 22898- 


027 


CS 31062- 


012 


CS 22880- 


-074 




[X/Fe]^ loge 




s.e. 


[X/Fe]^ loge 


s.e. 


[X/Fe]^ loge 


n}' 


s.e. 


[X/Fe]^ loge 




s.e. 


Fel 


-2.38 


5.12 


34 0.16 


-2.26 


5.24 


36 0.11 


-2.55 


4.95 


55 0.11 


-1.93 


5.57 


58 0.11 


Fe II 


-2.38 


5.12 


6 


0.20 


-2.25 


5.25 


4 0.16 


-2.55 


4.95 


3 


0.15 


-1.93 


5.57 


11 0.14 


Sr II 


+0.88 


1.42 


2 


0.43 


+0.92 


1.59 


2 0.41 


+0.30 


0.67 


2 


0.40 


+0.39 


1.38 


2 


0.40 


■\7" TT 

Y II 










+0.73 


0.71 


2 0.20 


+0.59 


0.27 


1 


0.20 


+0.16 


0.46 


2 


0.18 


7t- TT 
Zii 11 


+1.11 


1.34 


1 


0.23 


+1.01 


1.37 


2 0.17 


















Ba II 


+2.11 


1.95 


2 


0.17 


+2.23 


2.20 


3 0.14 


+1.98 


1.65 


1 


0.16 


+1.31 


1.60 


1 


0.16 


La II 


+1.69 


0.53 


6 


0.24 


+2.13 


1.10 


5 0.18 


+2.02 


0.69 


3 


0.18 


+1.07 


0.36 


3 


0.18 


Ce II 


+2.01: 


1.26: 


1 


0.22 


+2.13 


1.51 


18 0.13 


+2.12 


1.20 


7 


0.13 


+1.22 


0.92 


3 


0.14 


Nd II 


+2.01 


1.12 


2 


0.23 


+2.23 


1.47 


8 0.14 


+1.79 


0.73 


4 


0.14 


+1.20 


0.76 


1 


0.18 


Sm II 










+2.08 


0.82 


3 0.15 


















Eu II 


+1.73 


-0.10 


3 


0.20 


+1.88 


0.18 


3 0.15 


+1.62 


-0.38 


3 


0.14 


+0.5 


-0.8 


1 


0.17 


Dy II 










+1.78 


0.70 


2 0.17 










+0.60 


-0.16 


1 


0.19 


Er II 










+2.40 


1.12 


4 0.17 










+1.41 


0.45 


1 


0.19 


Pb II 


+3.3 


3.0 


2 


0.24 


+2.84 


2.65 


2 0.19 


+2.4 


1.9 


1 


0.19 


+1.9 


2.0 


1 


0.19 




CS 31062- 


050 


HD 196944 


CS 22942- 


019 


CS 30301- 


-015 




[X/Fe]^ 


loge 




s.c. 


[X/Fc]^ 


loge 


s.c. 


[X/Fe]^ 


loge 




s.e. 


[X/Fe]^ 


loge 




s.e. 


Fe I 


-2.31 


5.19 


26 0.14 


-2.25 


5.25 


76 0.19 


-2.64 


4.86 


18 0.14 


-2.64 


4.86 


30 0.18 


Fe II 


-2.33 


5.17 


3 


0.16 


-2.26 


5.24 


10 0.20 


-2.64 


4.86 


3 


0.15 


-2.63 


4.87 


9 


0.15 


Sr II 


+0.91 


1.51 


2 


0.18 


+0.84 


1.51 


3 0.22 


+1.7: 


2.0: 


2 


0.17 


+0.3: 


0.6: 


1 


0.23 


YII 










+0.56 


0.54 


7 0.22 


+1.58 


1.17 


2 


0.18 


+0.29 


-0.12 


4 


0.20 


Zr II 


+1.02 


1.31 


2 


0.16 


+0.66 


1.02 


6 0.20 


+1.69 


1.66 


3 


0.15 










Ball 
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Fig. 1. — Comparison of the observed spectra (dots) and synthetic spectra (hnes) near the Pb I 
4057.8A hne for HD 196944 and CS 31062-050. Sohd hnes show the synthetic spectra assuming 
the solar system Pb isotope ratio (see text). The Pb abundances assumed in the calculations are 
[Pb/Fe] = 1.6, 1.9, and 2.2 for HD 196944, and [Pb/Fe] = 2.6, 2.9, and 3.2 for CS 31062-050. The 
dotted lines show the synthetic spectra calculated with a single line approximation for [Pb/Fe] = 
1.9 and 3.1 for HD 196944 and CS 31062-050, respectively. 



Fig. 2. — The abundances of heavy elements as a function of atomic species for HD 196944 (upper) 
and CS 31062 050 (lower). The thick solid line indicates the main s-process component determined 
by Arlandini et al. (1999), while the thin line indicates the r-process component. The dotted line 
represents the total solar abundance adopted in Arlandini et al. (1999). All abundance patterns 
are normalized to the observed Ba abundances. 



Fig. 3. — [Pb/Ba] as a function of [Fe/H]. Filled circles show the results of the present work and 
of Aoki et al. (2002b) for LP 625 44 based on the Subaru/HDS spectra. The open circle shows 
the ratio of LP 706-7 determined by Norris, Ryan, &; Beers (1997) and Aoki et al. (2001) based 
on the AAT/UCLES spectrum. 
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